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INTRODUCTION 

Me  have  ahown  (l#21  that  orthoqonal  collocation  la  suitable 
not  only  for  the  solution  of  linear  diffusion  problew  but  also 
for  sisulatlon  of  tuo-dinensional  aodels.  Itius.  It  is  possible 
to  slaulate  edqe  effects  on  electroanalytlcal  eaperlnents  usinq 
this  fast  and  flesible  nuBerical  technique. 

Cyclic  voltaneetry  is  one  of  the  BM>st  rn— only  used 
electroanalytlcal  Methods  and  it  ntqhc  be  of  considerable 
Interest  to  calculate  the  effects  caused  by  edqe  diffusion  on  the 
i/E>curves.  Lines  and  Parker  (S)  have  observed  an  increase  of 
aCp  Mhen  uslnq  aicroelectrodes  and  attributed  this  to  edqe 
effects,  a  conclusion  which  wss  eonfirtsed  by  theoretical 
calculations  by  Heinse  |41.  Spelser  and  Rieker  (SI  found 
draoatlc  chanqes  in  cyclic  voltasuaoqrsMs  of  2*S-Di-^rt  butyl* 
4(4*diMethylaisinophenyU*phenol  I  in  unbuffered,  neutral 
acetonitrile  when  decreasinq  the  electrode  area  fros  0.29  cn^  to 
0.0022  ee^.  These  esperiaents  have  also  been  etplained  in  tents 
of  e^qe  diffusion,  up  to  now,  however,  it  was  not  possible  to 
aJMlajl*  the  coNpiet  aechaniai  descrlbinq  the  esldation  of  I 
under  the  conditions  of  a  finite  circular  disc  electrode. 


ihis  paper  expands  the  theory  of  two-diaensional  orthoqonal 
collocation  to  cyclic  voltawMtry  and  eechanians  aa  coaplex  aa 
the  oxidation  of  I.  It  will  only  deal  with  the  aetheisatics  of 
the  Method,  while  a  cotplete  nuterlcal  analysis  showinq  the 
usefulness  of  the  derived  equations  will  be  qlven  in  part  4. 


This  (Mrtition  has  b«en  us«d  because  the  analysis  (which  has  to 
be  done  to  validate  the  new  two-diae ns  tonal  collocation  ■ethodl 
requires  the  discussion  ot  aany  paraiMters.  Ue  point  out« 
however,  that  this  has  to  be  done  only  once,  while  extension  to 
•ore  cooplicated  aechanisas  is  obvious.  Also,  research  on  a 
coaparison  of  the  experiaental  results  with  1  and  nuserical 
siaulations  is  under  way. 

INITIAL  AND  BOUNOARV  CONDITIONS 

The  initial  and  boundary  conditions  have  been  qiven  by 
Nicholson  and  Shain  (A)  and  used  by  Speiser  and  Rieker  (7)  in  the 
contest  of  orthogonal  collocation  for  linear  diffusion,  pons  (8) 
siaulated  cyclic  voltaenograas  under  spherical  diffusion 
conditions  using  orthogonal  collocation. 

For  sinulation  of  edge  eftects  we  eapluy  our  two^iaensional 
•odel  of  a  finitSi  circular  electrode  of  radius  r^  labedded  into 
an  insulator  body  (1),  The  radius  of  the  insulator  aay  go  to 
infinity. 

He  assuwe  an  ec  •echanisAj 


Q  — »  C  A.  ‘ 


where  the  diffusion  controlled  electron  tcatisfer  is  described  by 
the  toreal  potential,  C*,  and  the  nunber  of  transferred 
electrons,  n.  The  chenical  step  (2)  proceeds  with  a  rate 
constant  k. 

Pick's  second  law  tor  the  substrate  A  takes  tl«e  torn 


where  t  is  t iau  and  s  and  r  are  the  distance  coordinate 
perpendicular  to  the  electrode  and  the  radial  coordinate 
respectively.  The  total  diffusion  is  espressed  as 
superiaposition  of  a  linear  end  a  cylindrical  tera. 

Por  species  B  end  C  aquations  14)  and  IS)  lioldi 

it  L  T  Jr  Jr*  J  ^ 

ra*c,  4  Jt,  1 

3t  ®  I  T  Jr  9r*  J  * 

Por  convenience  t)ie  diffueion  coefficient  of  ell  epeciea  ie  set 


Ite  assuee  that  A  ia  the  only  epeciee  present  in  the  solution 
at  the  beginning  of  the  esperiaent*  its  concentration  being  cj^* 
Thus,  the  initial  conditions  are 

t‘o,  >ts€>, 

It) 


(^()^,r,p)xt^(l,^,p)s.O  ' 

At  a  disiancv  far  froi  the  electrode  no  concentration  chenqea 
take  place  during  the  esperiaent.  Thie  yields  the  boundary 
conditions  at  s  •  •  and  r  *  • 
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functions  hav«  to  fulfill  tho  conditions  (20)  and  (21)  at  r  «  tq 


for  all  siwcies. 

t>o,  r»r,: 
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TRANSFORMATION  TO  A  DIMeNSlONLRSS  PROBLEM 


Squat  .ona  (1)»(22)  are  tranaforiMd  into  a  disenaionleas 
problea  by 
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(28) 


P.*r/r,  (^oirsr, 

R»(r- *■,)/(«-%)  (x-i-vn  (29b> 

(30) 

o<eAW 


Here  X  la  a  diaenalonleaa  distance  perpendicular  to  the 
electrode,  L  is  a  distance  froa  the  electrode  defined  such  that 
diffusion  effects  are  negligible  at  all  tiaes  tlie  c*  are 
diaensionless  concentrations,  R  and  R*  are  diswnsionlees 
coordinates  in  r-direetion  for  the  inner  and  outer  functions 
respectively,  uhile  M  is  defined  as  L  eseept  it  is  a  radial 
distance  froa  the  center  of  the  electrode,  llte  syebole  T*  end  e 
denote  t))e  dioensionless  sinuletion  tiae  and  rate  constant 
respectively. 

Thus,  ue  get  the  diMhsionless  differential  eguationst 
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AS  in  tbs  csss  of  tbs  cbroftosiipsrosstrlc  siaulstlons  Cl]  ^  bsv« 


tHD  dlssnslonlsst  psrsasfcers*  hsrs  $  snd  s".  9bs  psrssstsr  S*  is 
qU»n  by  espsrlsentsl  pstsMtsrs  snd  thus  ths  sUalstlons  ulll 
depend  on  0*  ss  e.q.  on  •. 

Ths  trsnsfoncd  initlsl  snd  boundsry  .onditlons  srs  qUen 
by 
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where  N,  end  il^  ere  the  nwebere  of  interior  eollocetion  points 
(chosen  ss  the  roots  of  an  orthoqonal  polynoolal  19))  In  the  i 
and  r  direction  respectively  and  and  are 

elements  of  natrices.  Hie  nunerical  values  of  those  eienents 
depend  only  on  the  decree  and  liind  of  the  polynonial  chosen  for 
collocation  (10).  thua,  if  the  ostrli  eieoents  have  once  been 
calculeted  they  nay  be  used  for  any  problem.  The  decree  of  the 
potynoailal  used  depends  on  the  desired  seeurscy  of  the 
simulation. 

Usln9  equations  (42)-<69)  ws  msy  discretite  the  psrtial 
dlfferentisl  equations  (!])•( )7)  in  teros  of  the  concentrations 
at  the  collocetlon  points,  ft  we  Introduce  the  boundary 
conditions  (43)«(40)  ws  srrlve  at  equations  (64)-(7U. 
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priMd  Mtrii  «l«MnU  correspond  to  the  outer 
diecretitetion  function.  As  in  the  cese  of  chronuesperoaetry  im 
hove  to  substitute  iaplicit  equations  for  the  unknown 
concentrations  at  the  boundaries.  The  calculations  of  the 
cl(X|,0«TM.  cJiXt.l.TM.  cJiXi.O.t')  Ci(Q,Rk,T'I  end 
^2t0*Rt«T*)  is  essentially  the  sane  as  for  chronoaeperonetric 
siaulations  and  we  refer  to  part  1  of  this  series  111  for  the 
derivations. 


The^  are  |  11  : 


ixifip  • '  i:  ( 

"  t‘i  «-i  *  ■ 


♦  C  c  I 


CH) 


(XiXj'J  * 

‘iU;  p,T') ~  ♦ 

<Oci^,r„c:,(jc.oj'h  f  T-) 

»c;,(x„ar7.  ^f^n,  V’J 

^*e.>  (HI 


4  I-  ,  V’ 

e  i  21  He  fci^*»^«  J  M 

^  V  #.»  f*J  *  -/ 


^  “•“* 

fn 


A  [ 

^•,(o.V’)--  -  7-  I  .:.(y,,Pi,T')] 

V*  I 


^f«4 


c|.(o.*',.r)  r  -  ^ •  c;,  (Xj.  f; J7 


C(.,,V,V).-  -  r  <,(x(,«vr) 


(03) 


%4  4:x 


(84) 


^f.A  “  Cf, 


(88) 


The  values  of  c^j(OiR)(,T')  and  Cri(0  »R|(,T* )  «ay  be  calculated 
froe  equations  (48)  and  (SO),  yleldlnq 
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How  we  a«y  substitute  (74)«(87)  Into  equations  (88)«(7l)  which 
qlves  us  finally  the  dlscretlsed  equations  deserlblnq  the  tine 
dependence  of  the  concentr 
at  the  collocation  points. 
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Pro*  (6$)  U  (ollowB  using  (74)  (77)  snd  (84) 
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stiiilarly  (68)  gives  with  (76)*  (  79)  snJ  (80) 

A  (Li  fi  ;y‘ 
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e*t  j 

Per  the  outer  function  describing  the  concentrstion 
species  A  it  holds  fros  (69)  •  (77)*  end  (81) 

♦  l-f VA*|ki  Hf  <*  Oshp 


If  MS  cosbins  (67)  with  (7S)«  (78)  end  (87).  we  get 
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Thuft.  the  syst««t  of  6  second  order  pertlel  differential 
equations  (32)'‘(3'7)  has  been  converted  into  a  set  of  3N^(Np  * 

Np>)  sieultaneous  ordinary  differential  equations*  (18)  and  (9U- 
(9S).  After  Inteqratinq  this  set  ue  qet  the  concentrations  of 
species  k,  B  and  C  at  the  collocation  points  (XpP|()  and  <X|« 
as  a  function  of  the  sinulation  tine  T' . 

One  of  the  distinctive  features  of  the  orthogonal 
collocation  approach  can  be  seen  fron  a  conparlson  of  equations 
(98)  and  I91}-(9S)«  Those  fornulae  are  very  sinllar  fron  a 
structural  point  of  vies,  rurtheroore*  coapared  to  the 
expressions  deacrlblnq  the  aodel  of  a  chronoaaperoaetrlc 
eiperiaent  (1| ,  there  are  only  few  changes  in  the  inner 
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functions,  while  the  outer  functions  ate  the  sane  both  for  the 
ixodels  of  cyclic  voltaxwietry  and  chronoaaperoMtry. 

CALCULATICM  OP  THE  CURRCMT 

Also,  the  calculation  of  the  current  Is  fully  analogous  to 
that  performed  for  the  chronoanperoaetrlc  aodel.  Again  we  have 
to  integrate  the  flux  over  the  part  between  0  and  on  the  r> 
axis 


O 


Transforming  the  right  hand  side  of  equation  (98)  Into  a 
dimensionless  form  we  get 

or,  using  Nicholson  and  Shaln*s  current  function  (8) 


The  integration  may  be  performed  exactly  by  a  quadrature  method 
(HI*  for  We  )(now  the  derivative  (  collocation 

points 
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where  the  el«»enta  of  the  vector  0  depend  only  on  thu  de^ri.*e 
kind  of  the  triel  function  in  r-direction. 


MODEL  MECHANISM  POA  OXIDATION  OP  t 


It  hes  been  shown  I12,J3|  that  the  oxidation  of  1  (denoted 
as  A)  say  be  described  by  an  ECE-type  aechanisn.  where  protons 
released  after  the  oxidation  step  are  picked  up  by  unoxidixed 
phenol  aolecuies  qiuihg  the  protonated  font  of  I,  HA*.  The 
latter  is  oxidised  at  more  positive  potentials  releasing  two 
protons  during  reaction.  The  end  product  of  both  reactions  is 
the  stable  phcnoxeniuA  cation.  E,  which  is  leost  probably  forned 
(ruN  the  neutral  phenosy  radical.  C,  in  a  hosogeneous  redox 
reaction  between  C  and  the  prisary  radical  cation  B  U-JJ  •  The 
intersediate  o  was  not  oxidised  further.  The  observation  of 
Changes  in  voitaMograss  of  1  with  electrode  size  has  been 
attributed  to  the  protonation  reaction  A-»MA*  which  causes  a 
phenol-ftee  layer  in  front  of  the  electrode  i 5| .  The  tine  after 
wtiich  phenol  «gain  reaches  the  electrode  surface  is  decreased  if 
•ore  substrate  diffuses  into  this  layei  dueto  edge  effects. 
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WiUt  Our  two>diaensional  sihulation  nodei  it  is  now  possible  to 
include  the  effects  of  a  finite  disc  into  the  calculations.  This 
•ay  also  serve  as  a  further  exaeple  for  the  application  of 
orthogonal  collocation  to  electroanaiyticai  problens.  In  this 
paper  we  show  the  derivation  of  the  expressions  for  dc*/dT*  for 
the  species  involved,  vhile  the  conparison  with  experieental 
results  will  be  reported  in  part  4  of  this  series. 

The  kinetic  diffusion  equations  for  the  above  node! 

•echanm  in  dieensionless  fom  are 


It* 


J'c* 

?K* 


E. 


“«*■«%  *».  “ii  ‘e* 


1102) 


^  ./1  ♦  4.’  ?:lj  c*  ^ ,» .... .. 


( 104) 


^  /l’ Jet., 

>T-  lx*  ■'r.  'Ip 


(  lOS) 


t 

1 


s  • 


/ 


u 


Jx*  R 


2^!; 

J«' 


<  106) 


<  107) 


9r  -r -jF *  R  IS  av  "-C<^ 

?r-  f  Jx«  a*r-«,-V;.-«v;.V.-*.‘.V:, 

I 

iiK.  fsS  ..  9V, 

»r  T  JX* 


(lOS) 


(  109) 


UlO) 


2^, A  1^’^  5'ii  JcJ 


9j% 

»r 


^  9'»»f  j<  afi5»» _ .  .  . ,« 


ilis.A  »•  V* .  .. 

JT'  '  9X*  »R’  »(r*  V 

’'ft  A-’^:, 

>T-  »x*  9V*^ 


( 114) 


!!e..a’j[A  7tf,  9\f,  ^  ^ 

TV  -1^  9x‘  {]ir*ljPfi'  9r  ■'''  ST*  * 


wher* 

vV*'  '  '*'*> 

«i'‘4»«aA  (  1171 


and 


.r.«-a.,t;/o.  (M»k; 


1  • 


I 


In  addition  to  the  initial  conditions  (41)  and  (42)  we 
assujse 
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The  boundary  conditions  include  now 
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CONCLUSION 

We  showed  in  this  paper  the  use  ot  two-dinensionsl 
collocation  to  spatially  discretite  partial  difCerentisl 
equations  describinq  a  oodel  of  cyclic  voltanactric  experlisents 
under  the  influence  of  edqe  effects.  Hie  forisal  slsiilarity  of 
descretited  equations  for  even  very  conplicated  nechanisns  has 
been  denonstrated  by  conslderinq  the  oxidation  echeise  of  phenol 
I.  This  slsiilarity  sinplifies  the  wrltinq  and  testlnq  of 
conputer  proqrans  used  to  calculate  niMerical  sisulation 
results.  In  fact*  e  ehanqe  in  siechanian  nornally  only  requires 
the  eschanyo  of  a  few  PONTflMi  stateoents  while  the  nain  body  of 
the  proqrao  is  untouched.  Me  see  this  as  one  of  the  qreat 
advantaqes  of  our  Method  which  outweighs  the  difficulties 
encountered  while  assessing  the  features  and  behavior  of  a  new 
type  of  a  Model*  as  this  two^dinenaional  one.  A  thorough 
discussion  of  the  nuMerical  results  for  cyclic  voltaMMetry 
sinulatlons  is  given  In  pert  4. 
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